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Abstract 

The gauge-extended U{l)c x SU{2)l x U{l)ij^ x U{1)l model has the attractive property of 
elevating the two major global symmetries of the standard model (baryon number B and lepton 
number L) to local gauge symmetries. The U{1)l symmetry prevents the generation of Majorana 
masses, leading to three superweakly interacting right-handed neutrinos. This also renders a B — L 
symmetry non-anomalous. We show that the supcrweak interactions of these Dirac states (through 
their coupling to the TeV-scale B — L gauge boson) permit right-handed neutrino decoupling just 
above the QCD phase transition: 175 MeV < T^^"" < 250 MeV. In this transitional region, the 
residual temperature ratio between and vr generates extra relativistic degrees of freedom at 
BBN and at the CMB epochs. Consistency (within la) with both WMAP 7-year data and the 
most recent estimate of the primordial ^He mass fraction is achieved for 3 TeV < Mb-l < 6 TeV. 
The model is fully predictive, and can be confronted with dijet and dilepton data (or lack thereof) 
from LHC7 and, eventually, LHC14. 



Heavy neutral vector gauge bosons (Z'^s) are ubiquitous in extensions of the standard 
model (SM) [T], often including a gauged B — L symmetry. This U{1) symmetry is non- 
anomalous if the three left-handed Weyl neutrinos are accompanied by three right-handed 
neutrinos. The Z' masses are a priori open parameters - not determined by the low energy 
effective theory - but subject to recent experimental bounds (M^/ > 3 TeV) from searches 
of dilepton |2] and dijet |3] events in the 7 TeV run of the Large Hadron Collider (LHC7). 
In this Letter we re-examine some critical cosmological issues surrounding the presence of 
the six additional neutrino degrees of freedom [3] correlated to the presence of a Z' in 
our dynamical model which is coupled to B — L. These considerations, when viewed in 
the context of most recent data collected by the Wilkinson Microwave Anisotropy Probe 
(WMAP) [5j, are found to constrain the mass of the Z' to an interestingly narrow band, 
which will be directly probed by LHC14. 

For a good part of the past two decades, big-bang nucleosynthesis (BBN) provided the 
best inference of the radiation content of the universe. The time-dependent quantity being 
the neutron abundance at t > r„, which regulates the primordial fraction of baryonic mass 
in ^He, 
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where 5Nl^ = N^^ — 3 is the effective number of extra (non-SM) light neutrino species, 
r„ is the neutron half life, and rj is the ratio of the baryon number density to the photon 
number density [6] . The observationally- inferred primordial fractions of baryonic mass in ^He 
(Fp = 0.2472±0.0012, Yp = 0.2516±0.0011, Yp = 0.2477±0.0029, and Fp = 0.240±0.006 [7J) 
have been constantly favoring A^,^^ ^ 3 [8j. Out of the blue, two recent independent studies 
determined Yp = 0.2565 ± O.OOl(stat) ± 0.005(syst) and Yp = 0.2561 ± 0.011 [9j. For 
r„ = 885.4 ± 0.9 s and r„ = 878.5 ± 0.8 s, the updated effective number of light neutrino 
species is reported as A^^^ = S.GStofo (2a) and A^^^ = 3.80l[!;?[! {2a), respectively. 

Very recently, in support of these trends, observations of the cosmic microwave back- 
ground (CMB) anisotropics and the large-scale structure distribution have allowed a probe 
of A^^*^ at the CMB decoupling epoch with unprecedented precision. The relativistic par- 
ticles that stream freely influence the CMB in two ways: (1) their energy density alters 
the matter-radiation equality epoch, and (2) their anisotropic stress acts as an additional 
source for the gravitational potential via Einstein's equations. Hence, the number of light 
relativistic species becomes a function of the matter density {Q^h^) and the redshift of 
matter-radiation equality (^cq), 
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where fl^h"^ = 2.469 x 10^^ is the present-day photon energy density (for Tcmb = 2.725 K) 
and the scaled Hubble parameter h is defined hy H = 100 h km Mpc~^ [lOj. The variation 
in A^,*;^ reads 
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The equality redshift is one of the fundamental observables that one can extract from the 
CMB power spectrum. More specifically, WMAP data constrain Zeq mainly from the height 
of the third acoustic peak relative to the first peak |5j. The fractional error in Qmh"^ is 
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determined using external data: the latest distance measurements from the Baryon Acoustic 
Oscillations (BAO) in the distribution of galaxies [11] and precise measurements of the 
Hubble constant Hq [T2]. The parameter constraints from the combination of WMAP 7- 
year data, BAO, and lead to Nf = 4.34+°|^ (68%CL) [S'j. 

In summary, though uncertainties remain large, the most recent cosmological observations 
show a consistent preference for additional relativistic degrees of freedom (r.d.o.f.) during 
BBN and the CMB epochs. We take these hints as motivation for our analysis, which 
consists of the following tasks: (1) to present a model in which the additional r.d.o.f. are 
three flavors of light right-handed neutrinos which interact with the SM fermions via the 
exchange of heavy vector fields; (2) to suppress the six additional fermionic r.d.o.f. to levels 
in compliance with BBN and CMB. This is accomplished by imposing the decoupling of i/^'s 
from the plasma early enough so that they undergo incomplete reheating during the QCD 
phase transition; and late enough so as to leave an excess neutrino density suggested by 
the data [13]. These requirements strongly constrain the masses of the heavy vector fields. 
Together with the couplings, which are determined in accord with other considerations, the 
model is fully predictive, and can be confronted with dijet and dilepton data (or lack thereof) 
from LHC7 and, eventually, LHC14. 

An economic choice of the model to implement the task outlined above is based on 
the gauge-extended sector U{3)c x SU{2)l x f/(l)/^ x U{1)l [H]. The resulting f/(l) 
content gauges the baryon number B [with U{1)b C U{3)c], the lepton number L, and a 
third additional abelian charge J/j which acts as the third isospin component of an SU{2)fi. 
The usual electroweak hypercharge is a linear combination of these three U{1) charges: 
Y = ^{B — L) + Ir. The matter fields consist of six sets (labeled by an index i = 1 — 6) 
of Weyl fermion-antifermion pairs: {Ur, Dr, Ll, Er, Ql, Nr). The field Nr_ is the right- 
handed neutrino (and left-handed antineutrino) accompanying the fields in the set Er, with 
mass ~ 1 eV. The gauging of lepton number precludes the presence of a seesaw for generating 
Majorana neutrino masses. In addition to the SM interactions, these fields experience two 
f/(l) gauge interactions mediated by two associated vector bosons {Z' and Z") whose masses 
lie well above a TeV. 

The initially free parameters consist of three couplings gB,gL,giR- These are augmented 
by three Euler angles to allow for a field rotation to coupling diagonal in hypercharge. This 
diagonalization fixes two of the angles and the orthogonal nature of the rotation introduces 
one constraint on the couplings P{gY, gs, gL, gi^) = 0. The baryon number coupling gs is 
fixed to be \/3/2 of the non-abelian SU{3) coupling at the scale of U{N) unification, and 
is therefore determined at all energies through RG running. This leaves one free angle and 
two couplings with one constraint. The two remaining degrees of freedom allow a further 
rotation leaving Z' to couple to B at 90% and Z" to couple to 5 — L at 99%.^ The U{1) 
quantum numbers and the physical couplings of the Z' and Z" to the fermion fields are given 
in Table |lj These couplings, which have been computed elsewhere [H], are functions of the 
charge assignments, the U{1) gauge couplings, and the mixing angles. All fields in a given 
set have a common g', g" couplings.^ 

The model as described enjoys distinct advantages: (1) Gauging of the anomalous B 
and its cancellation by generalized Green-Schwarz mechanism (which leaves i? as a global 



^ Although not generally appreciated, it is important to note that a 100% coupling of the Z' and Z" to B 

and B — respectively, is possible only if the U{1) gauge coupling constants are equal. 
^ Our couplings are consistent with the bounds presented in |15j from a variety of experimental constraints. 
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TABLE I: Quantum numbers of chiral fermions and their couplings to Z' and Z" gauge bosons. 



Name 


Representation 


B 


L 


Ir 


Y 


g' 


g" 


Vr 


(3,1) 


1 

3 





1 

2 


2 
3 


0.368 


-0.028 


Dr 


(3,1) 


1 
3 





1 
2 


1 

3 


0.368 


-0.209 


Ll 


(1,2) 





1 





1 

2 


0.143 


0.143 


Er 


(1,1) 





1 


1 

2 


-1 


0.142 


0.262 


Ql 


(3,2) 


1 

3 








1 

6 


0.368 


-0.119 


Nr 


(1,1) 





1 


1 

2 





0.143 


0.443 



symmetry) prevents proton decay. (2) The presence of Nr renders B — L non- anomalous. 
This has been appeahng for minimal extension of SM at the TeV-scale. For example, the 
mass growth of Z" can occur via a conventional Higgs mechanism at TeV without relying 
on possible Planck scale physics. (3) By inspection of Table |l] the charges B, L, and Ir 
are mutually orthogonal in the fermion space. This will maintain the othogonality relation 
P = to one loop without inducing kinetic mixing [13] . 

We begin by first establishing, in a model independent manner, the range of decoupling 
temperatures implied by the the BBN and CMB analyses. For this work, the physics of 
interest will be taking place at energies in the region of the QCD phase transition, so that 
we will restrict ourselves to the following fermionic fields, and their contribution to r.d.o.f.: 



[3ur] + [3dR] + [3ul + eL + Hl] + [cr + /iij] + [Sul + 3^^] + [3ur] . 



(4) 



This amounts to 22 Weyl fields, translating to 44 fermionic r.d.o.f. 

Next, in line with our stated plan, we use the data estimate to calculate the range of 
decoupling temperature. The effective number of neutrino species contributing to r.d.o.f. 
can be written as 
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therefore, taking into account the isentropic heating of the rest of the plasma between vr 
decoupling temperature Tjec and the end of the reheating phase. 



eff 



U(Tc 



doc J 



4/3 



(6) 



where Tend is the temperature at the end of the reheating phase, and A^(T) = r(T) (A^b + |A^f) 
is the effective number of r.d.o.f. at temperature T, with A^b = 2 for each real vector field 
and A^F = 2 for each spin-^ Weyl field. The coefficient r(T) is unity for the lepton and 
photon contributions, and is the ratio s(T)/ssb for the quark-gluon plasma. Here s(T)(ssb) 
is the actual (ideal Stefan-Bolzmann) entropy. Hence A^(Tdec) = 37 r(Tdec) + 14.25. We take 
A^(Tend) = 10.75 reflecting (e^ + + + elv^L + ^eR + ^^/.l + v^lR + VtL + VtR + 7l + Ir)- 
We consistently omit vr in considering the thermodynamics part of the discussion, but will 
include it when dealing with expansion. As stated in the introduction 
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so the excess r.d.o.f. will lie within la of the central value of each experiment if 0.46 < 
6Nf < 1.08. From Eqs. Q and (Q, the allowable range for is 

23 < iV(Tdec) < 44 . (8) 

This is achieved for 0.24 < r(Tdec) < 0.80. By comparing to Fig. 8 in Ref. [16], this can be 
translated into temperature range 

175 MeV < Tdec < 250 MeV , (9) 

with the lower temperature coinciding with the region of most rapid rise of the entropy. 
Thus, the data implies that the decoupling takes place during QCD phase transition. 

We now turn to use our model in conjunction with the decoupling condition to constrain 
its parameters. To this end we calculate the interaction rate r(T) for a right-handed neutrino 
and determine Tdec from the plasma via the prescription 

r(rdec) = i/(Tdec) . (10) 

Let fl be a single species of Weyl fermion, representing the two r.d.o.f. {fl, /^}, where the 
superscript indicates bins i = 3,5. Similarly /|j G /!}, for i = 1,2,4,6. Notice that 
the subscripts L, R denote the actual helicities of the massless particles in question, not the 
chirality of the fields. With this said, we may write the amplitude for fl scattering 

^ {Mpi)rLiP2) ^ MP3)fLiP4)) = ^[uiP3)ra + 7Mpi)MP4h,ii-i>iP2)] • (h) 

The other 3 amplitudes are obtained by the crossing substitutions in the second square 
bracket; for scattering from 

fk ^ w(P4)7m(1 + 75)w(P2) (12) 

The cross sections for the four scattering processes (no average over helicities) are 

a [vnfi ^ vnJl) = {unfn ^ '^Rfk) = ~ (for bins ^ = 3, 5) (13) 

and 

^ ^ = IcrMk -> lyRfk) = (forbinsz = 1,2,4,6) . (14) 

In addition to these scattering processes, the interacts with the plasma through the 
annihilation processes: i^/jz/^ — )■ for bins i = 3,5, and I'Ri^L fRpLi bins i = 

1, 2, 4, 6. These all yield cross sections 2Gfs/ (37r) due to forward and backward suppression. 
Assuming all chemical potentials to be zero, the plasma will have an equal number density 
n(T) = 0.0913T^, for each fermion r.d.o.f. Thus, 

r^^\T) = n{T) (y2a,is)vMM\ , (15) 
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where um = 1 — cos6'i2 is the MoUer velocity, s = 2kik2{l — 005^12) is the square of the 
center-of-mass energy, and TVi is the multiphcty of Weyl fields in each bin (e.g., for i = 
3, As = 3 + 2 = 5). The scattering cross section is given by 

- 4 2C'^ s 

= ^Ml ^ ^Rfl) + ^ ^Rfk) = foreachz = 1, . . . , 6 ; (16) 

O 71 



similarly. 



1 2(^2 s 

aris) = a{uRUL ^ flfn + 1^1) = for each z = 1, . . . 6 . (17) 



Since s = 2kik2{l — cosOu) and vm = 1 — cos ^12, we perform an approximate angular 
average ((1 — cos ^12)^) = 4/3, followed by a thermal averaging {2kik2) = 2(3.15^ T^) to give 



Tscat 



(T) = ^2 (3.15T)2(0.0919T3) fj^G^WJ ~2.05G'^gT^ (18) 



From (16), (17), and (18) 



pann^y^ = ir'^'=^*(T) ~ 0.50 G^g . (19) 

Each of the Gi is given by the sum of the contributions from Z' and Z" exchange, 

A _ 96 9i I 96 9i ,c^n.\ 
V2 " Ml, + M|„ • ^^""^ 

The Hubble expansion parameter during this time is 

H{T) = 1.66 (iV(T))^/2 TVMpi , (21) 

where Mpi is the Planck mass. Since the quark-gluon energy density in the plasma has a 
similar T dependence to that of the entropy (see Fig. 7 in [16]), we take N{T) = 37 r(T) + 
19.5, so that H{T) = 0.82 x 12.5 T^/Mpi. (The first factor provides an average for r(T) 
over the temperature region, and we have now included the six upt r.d.o.f.) Since F oc 
and H ~ T^, it is clear that if at some temperature Tdec, -f^(7dec) = rj(^dec)5 the ratio 



T/H will fall rapidly on further cooling. Thus from (10) and (21) the equation determining 
Tdec depends on: (1) whether we need to preserve the absence of a chemical potential, or 
(2) whether we need simply to mantain physical equilibrium. The decoupling condition in 
these two cases is: (1) F-'^(Tdec) = i^(Tdec) and (2) F-^t(Tdec) + F-"(Tdec) = i^(Tdec); or 
numerically: (1) 

0.50 Gls T|ec = 10.4 T^ec/Mpi Tl, = 20.8 (G.^Mpi)"^ , (22) 

and (2) 

2.50 Tl^ = 10.4 TlJMp, ^ T],, = 4.1 (Ge^^Pi)"' ■ (23) 

Tdec as determined from these equations must lie in the band (|9]). 

Since all freedom of determining coupling constant and mixing angles has been exercised, 
there remains only constraints on the possible values of Mz' and Mz"- Our results are 
encapsulated in Fig. [T| and along with other aspects of this work are summarized in these 
concluding remarks: 
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• In this Letter, we develop a dynamic explanation of recent hints that the relativistic 
component of the energy during the era of last scattering is equivalent to about 1 extra 
Weyl neutrino. 

• We work within the context of a specific (string based) model with 3 U{1) gauge 
symmetries, originally coupled to baryon number B, lepton number L , and a 3rd 
component of right-handed isospin I^. We find that rotation of the gauge fields to a 
basis exactly diagonal in hypercharge Y, and very nearly diagonal in B — L and B 
fixes all the mixing angles and the gauge couplings. Of course, of most significance for 
this work, requiring that the B — L current be anomaly free, implies the existence of 
3 right-handed Weyl neutrinos. 

• We then find that for certain ranges of Mb and Mb-l the decoupling of the vr^ occurs 
during the course of the QCD phase transition, just so that they are only partially 
reheated compared to the z/^^'s — the desired outcome. 

• To carry out this program, we needed to make use of some high statistics lattice 
simulations of a QCD plasma in the hot phase, especially the behavior of the entropy 
during the changeover. 

• Since our aim is to match the data, which has lower and upper bounds on the neutrino 
"excess", we obtain corresponding upper and lower bounds on the gauge field masses. 
Roughly speaking, if decoupling requires a freezout of the annihilation channel (loss 
of chemical equilibrium), then 3 TeV < Mb-l < 4 TeV. If thermal equilibrium via 
scattering is sufficient, then 4.5 TeV < Mb-l < 6 TeV. These are ranges to be probed 
at LHC14. 

• Finally, a remark about the model: the gauging of B allows a global conservation of 
baryon number. The gauging of L brings with it the loss of a heavy Majorana for 
the seesaw model, as well as for leptogenesis through the decay of this particle. Thus, 
along with all its companion fields, the neutrino it is a Dirac particle, with the small 
mass originating through small Yukawa. 

L.A.A. is supported by the U.S. National Science Foundation (NSF) under CAREER 
Award PHY-1053663. H.G. is supported by NSF Grant PHY-0757959. 
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FIG. 1: Dark shaded areas show region aUowed from decouphng requirements to accommodate 
WMAP and BBN data. Light shaded regions indicate the masses excluded by the LHC7 dijet 
searches. The dark shaded areas in the left and right figures pertain to chemical and thermal 
equilibrium, respectively. These two estimates should serve to bracket the size of the actual effect. 
The designation of B corresponds to Z' and B — L to Z" . 
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